Galilean invariance and homogeneous anisotropic randomly stirred flows 
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The Ward-Takahashi (WT) identities for incompressible flow implied by Galilean invariance are 
derived for the randomly forced Navier-Stokes equation (NSE), in which both the mean and fluc- 
tuating velocity components are explicitly present. The consequences of Galilean invariance for the 
vertex renormalization are drawn from this identity. 
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The transformation properties of the Navier-Stokes 
equation (NSE) play an important part in both the 
physics of turbulence as well as in practical aspects of 
turbulence modelling 0. Among the invariance proper- 
ties or symmetries of the NSE, Galilean invariance is per- 
haps the most frequently cited in theoretical approaches 
to turbulence 0, 0, 0, IE El ■ 

For fluctuating flows, a decomposition of the instanta- 
neous fluid velocity into a mean plus a component fluctu- 
ating relative to the mean (the Reynolds decomposition) 
shows that the mean and fluctuating velocities respond to 
a Galilean transformation in different ways. This would 
seem to be the most general physical significance one 
could deduce by studying the behavior of the NSE under 
Galilean transformations. On the other hand, because 
the Galilean transformation is a symmetry of the NSE, 
it leads to nontrivial and exact identities among various 
correlation functions. These so-called Ward-Takahashi 
(WT) identities make two types of important statements: 
(a) In quantum field theory (QFT), they relate formally 
divergent parts of different correlation functions. Ward- 
Takahashi identities reduce the number of counterterms 
or renormalization constants; and (b) in QFT as well as 
in field theoretic approaches to condensed matter prob- 
lems, they relate different physical quantities. For ex- 
ample, in hydrodynamics they relate the eddy viscosity 
and the numerical prefactor multiplying the nonlinear 
convective term (the vertex). Thus, the deeper phys- 
ical significance of the Galilean transformation is that 
it constrains the values of the physical parameters ap- 
pearing in the NSE, and in particular, it leads to the oft 
cited non-renormalization of the vertex of the NSE in the 
large-scale flow regime . 

Recently however it was claimed that Galilean invari- 
ance does not at all constrain the vertex renormalization 
, thus putting into serious question the widely cited key 
early papers in turbulence research 0,0- In 0, a distinc- 
tion is drawn between the NSE for the instantaneous ve- 
locity field V{ and the NSE for the fluctuating velocity Ui 
in the comoving frame of the mean constant velocity, for 



spatially homogeneous flow. These velocities are defined 
through the Reynolds decomposition of the velocity into 
its mean plus fluctuations: Vi(x.,t) = (Vj(x, t)) +Ui(x, t), 
where the angular brackets denote the ensemble aver- 
age pj. The essential observation made in is that 
Galilean invariance of the equation for the fluctuating 
velocity it, is established by transforming the mean ve- 
locity only. No previous field theoretic derivations of 
the vertex non-renormalization result have distinguished 
between the mean and fluctuating velocity components, 
since they treat the instantaneous velocity field at zero 
mean velocity, in other words, the case of a homogeneous 
and isotropic flow. In order to address specifically the 
concern in , in which the mean velocity plays a distin- 
guished role, one should examine carefully the derivation 
and physical consequence of the pertinent WT identities, 
keeping explicit the decomposition of the instantaneous 
velocity field into its fluctuating and generally nonzero 
mean components. In short, one needs to treat the case 
of a homogeneous but anisotropic fluid. The purpose 
of this Brief Report is to carry this out in some detail, 
paying careful attention to the rather different ways the 
mean and fluctuating velocities are affected by a Galilean 
transformation. The issue of vertex renormalization will 
be addressed within this framework. 

We start with the Reynolds decomposition to express 
the NSE in terms of the constant mean velocity and the 
fluctuation about the mean. So we take Vi — Ki + it,;, 
where Ki — (Vi) is the constant mean velocity and Ui 
is the fluctuation about the mean. By definition, (itj) = 
0. As demonstrated in 0, the equation of motion for a 
constant mean velocity Ki vanishes identically, resulting 
in the NSE for the fluctuating velocity u, . To this we add 
a stirring force fi and express the equation in solenoidal 
form: 



duj 
dt 



dui 



d{un 



dxi 



^oVV + U, (1) 



where the projection operator is Pij ( V) 
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Vj^yVj). fi is Gaussian with zero mean (fi) = and 
time and space translation invariant: (/i(x, t)fj(y, r)) = 
Dij (x— y)5(t — t) . The viscosity is written with subscript 
zero, i>o, to remind us that this is the bare (unrenormal- 
ized) molecular viscosity that appears in the fluid equa- 
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tion of motion. We provide the nonlinear term in Eq. 
l[TJl with a parameter Ao = 1. This is a well-known con- 
venient bookkeeping device used in perturbation theory, 
in which typically, solutions of Eq. are developed 
formally in powers of Ao- 

As we are interested in the correlation functions im- 
plied by Eq. iHTl, we make use of dynamic functional 
methods [8J, |9|, llCM. Ill| . The generating functional for the 
correlation functions associated with Eq. is given by 



Z[J,H] K = [[Du][Da]exp{-S[u,a] 



K 



with an analogous expansion for a. This implies that the 
functional Z in Eq. transforms as Z — » Z + 5Z, where 

SZ = <( / tJi(Sc-V)ui+m(Sc-V)a l -J-Sc s )) t = 0. (5) 

The angular brackets with bullet denote the average 
taken over the turbulent ensemble in the presence of the 
source terms. As Eq. JSJ holds for all Sc, it follows that Z 
satisfies the functional differential equation (the integral 
is over x, t) 



f if/ tJiV,J- + t^iViJ- - Jj)z[J, S] K = 0. 

+ / dxdi(J-u + J-K + E-<7)L (2) \J J SJi 3 5^i 3 ) L ' JA 



(6) 



where the action S is 



S[u,a] K = J dxdt(±J dya i {x,t)D ij (x-y)a j (y,t) 



/du k du k d{uiUj) 



dxi 



(3) 



The field a is conjugate to the velocity Ui, and is the ana- 
log of the noncommuting operator introduced by Martin, 
Siggia and Rose (MSR) in their operator formalism of 
classical statistical dynamics jl^. J and X are arbitrary 
space and time dependent source functions for V = u+K 
and a, respectively. It is important to note that both the 
mean and fluctuating velocities couple to the same exter- 
nal source. A straightforward and convincing way to see 
that this must be so is to start from the forced NSE for 
the instantaneous velocity V, and set up its correspond- 
ing action and dynamic functional. Couple Vi to a source 
term J^. Insert the Reynolds decomposition Vi — Ki + Ui 
into the measure, action and into the source terms. The 
incompressibility of the flow implies the fluctuations are 
themselves incompressible: ViVi = => V^Ui = 0. Then 
Eqs. and |JS} follow immediately from the dynamic 
functional and action for the instantaneous velocity. 

We now subject the generating functional Z to a 
Galilean transformation. Consider a second primed 
frame moving with constant velocity c with respect to 
the unprimed frame. Then the relations between time 
and coordinates of an event in the unprimed and primed 
frames are t = t' and x = ct + x' while the velocities 
and conjugate field transform as Ki — Ci + K[, Ui(x, t) — 
u'i(x!, t') and <7i(x, t) — er-(x', t'). Note that the mean and 
fluctuating velocities are affected in different ways. The 
conjugate field a has dimensions of acceleration and so 
transforms as indicated above. Using the time and coor- 
dinate transformations as well as these rules, it is easy to 
verify that S[u' , o~']k' = S[u, <j]k and the functional inte- 
gral measures are invariant: [Du'] = [Du], [Da'] = [Da]. 

To obtain the Ward-Takahashi identity, consider an 
infinitesimal Galilean transformation and Taylor expand 
the fluctuating velocity out to first order in the relative 
velocity 5c as follows: 

u(x + 5ct,t) = u(x,t) +t(Sc- V)u(x,t), (4) 



This makes use of the fact that the mean velocity Ki is 
constant in homogeneous turbulence, hence Vj(Ki) m — 
Vj--Ki(l). - KiVjil). - 0, and so V^^[J,S] K = 
V 3 -(ttj + Ki). = Vj(u/}. + Vj{Ki). = Vj{ui}.. This 
identity Eq. for the functional Z can be used to 
obtain exact relationships among the correlation func- 
tions associated with the random velocity field in its 
Reynolds decomposed form. Of course, Z generates 
both the disconnected and connected correlation func- 
tions. For the purposes of renormalization, it is custom- 
ary to work instead with the connected proper vertices, 
also known as the one-particle irreducible (1PI) Green 
functions 0] . It is important to point out that at lowest 
order, or tree-level, the connected proper vertex func- 
tions coincide with the vertices of the original (bare) ac- 
tion S ©■ To obtain the generating functional for 1PI 
vertices, we follow standard practice and introduce the 
generating functional of the connected Green functions 
W — InZ, and then carry out a Legendre transform [l3| 
r[(V)„ (a).] = -W[J, S] + / J-(V). + £»., where 



(Vk). 



5W 



SW 5T 
S(V k ). 



J k , 



6T 



S{ak}» 



(7) 

Note of course that (Vk), = (tife). + Kk- From the defi- 
nitions of W and T and using Eq. J7J, we obtain imme- 
diately the Ward-Takahashi identity satisfied by T: 



(tVi(Vi). 



ST 



-tVi(ai). 



ST 



ST 



) = 0. (8) 



The great utility of the WT identity Eq. © is that it pro- 
vides exact, non-perturbative equations relating various 
proper vertices (correlation functions) associated with 
the NSE equation. At lowest order (i.e., tree level), these 
proper vertices can be read off directly from the action S 
(see below). Here, and in conformity with the Reynolds 
decomposition, the mean value of the velocity field in the 
absence of source terms is a constant (and nonzero) vec- 
tor: (Vi(x,t)} = 5W[3, H]/6Ji | j= E =o = Ki. A nonzero 
expectation value for the fluid velocity in the limit of van- 
ishing external sources is reminiscent of symmetry break- 
ing in field theory. Indeed, the existence of a constant K j 
breaks the fluid isotropy since a preferred direction is be- 
ing singled out. In this case, the correct Taylor series 
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representation of the effective action T is given by 0] 

1 

ml 

ml,m2-l 
,i2, . . . ,iml ,jl,j2,. . . .jrn2 ( 



((V5mi(ml)>. -if„„i)(a J i(l))....(CT jm2 (m2)). (9) 



r[(V}.,(<7).]= f] — i- / 
mV.m2\ I 

™ i ™ o i « 



[dx rft] 



-^il,i2,... jml,.7'l,.?'2,...,.7'm2 (!))• -^il)-" X 



where 



i(ml,m2) 
il,i2,...,iml jl j'2,...jm2 



£ ml+m2 r[(v}.».] 



5(y i i(l)>....(5(Kmi(ml)).J((7 J -i(l))....5<o-j m 2(m2)). 



<V>.=K 
<<t).=0 



(10) 



We employ a condensed notation to indicate the de- 
pendence on space and time, thus, e.g., (1) = (xi,ii), 
(ml) = (x m i,t m i) and [chcdt] stands for the volume ele- 
ment for all ml + m2 pairs of space and time points to be 
integrated over in each term of the sum in Eq. ©. The 
factors Eq. I|10(l for each ml and m2, correspond to spe- 
cific proper vertex with ml factors of the velocity and ml 
factors of the conjugate field. So, for example, the inverse 
response function and proper vertex associated with the 
nonlinear convective term in the NSE are (suppressing 

the dependence on the space and time arguments) F^' 1 ^, 

and rj;V"\ respectively. 

The pertinent WT identity that we seek is obtained by 
inserting Eq. @ into Eq. JHJ differentiating the latter 
with respect to xlw , s ., u t? — -, s f „ u , using the defini- 

^ a<Vfc(y,t'))« <5<0j(w,t")). ' b 

tion in Eq. I|10|l and then setting (V). = K and (a), = 0. 
Doing so yields 

= J <txdtr$\x,t,y,if,w,ir), (ii) 

which follows after an integration by parts (we adopt 
vanishing boundary conditions) and using the fact that 
spatial homogeneity implies proper vertices are transla- 
tional invariant in configuration space. For our final step, 
we express this identity in terms of wavenumber (or, mo- 
mentum) and frequency space (k, oS), by means of the 
Fourier transform (FT). Recall that for stationary flow, 
the proper vertices can depend only on time differences. 
Thus, the FT of Eq. JTTJ yields 




(12) 

At zero-loop order (indicated via the zero superscript), 
and from Eq. © we can read off directly the 1PI func- 
tions that are related by this WT identity, namely 

o (M) / ,\ 

T tj = (-iLu + iKjk-j +v a k 2 )P lj (k) (13) 

(2,i) 

T ijk = *AofciPjfc(k). (14) 



It is important to realize that the ^ -dependence of the 
bare response function H13[) follows from applying the 
Reynolds decomposition from the outset to the NSE, 
which leads to JIJ and ©. This expression cannot be 
obtained from applying a Galilean transformation to the 
response function for an isotropic fluid (Kj — 0). Insert- 
ing Eqs. (|13I14|) into Eq. ljT2*j) implies that A = 1. In 
other words, Galilean invariancc requires that the book- 
keeping parameter introduced in Eq. Q must be iden- 
tically unity. This result is clearly independent of the 
mean velocity Ki. To go further, we make the reason- 
able assumption that some of the perturbative correc- 
tions to the solution of Eq. can be absorbed into re- 
definitions of the parameters appearing in the equation 
of motion. This is to say, we assume the NSE is partially 
form-invariant with respect to the renormalization aris- 
ing from the fluctuations inherent in the randomly stirred 
ensemble. This renormalization leads to frequency and 
wavevector dependent viscosity and a "mass" term S in 
the inverse response function. The vertex associated with 
the convective term in the NSE can also in principle 
be modified or suffer renormalization, leading to a fre- 
quency and wavevector dependent term which we denote 
A. The most general terms one can write down which 
maintain the same tensorial structure as their unrenor- 
malized counterparts Eq. \1'6\ and Eq. (|14|) are 

rf/^k^-k, -u) = 

( - iu + v{uj, k)k 2 + S(w, k)jP jk (k) (15) 

r^ 1} (0, 0; k, w; -k, -w) = ik t P jk (k) + A ijk (u, k), 

(16) 

where in accordance with Eq. (|12fl and Eqs. i|13ll4|) . we 
have set Ao = 1 in Eq. (fTB|> . The relation of Eq. ijTB|) 
and Eq. (|To)) to Eq. (fP3*)l and Eq. ijTil) , respectively, is as 
follows. Imagine setting the nonlinear term in the NSE 
to zero by formally taking Ao — > 0. Then v(u>, k) — > Uq, 
H(u>,k) — > iKjkj and Aij k (u>, k) — * 0. That is, these 
terms appear as a result of the perturbation expansion, 
which is developed in powers of Ao. Finally, inserting 
these expressions Eqs. <|15I16[) into the Ward-Takahashi 
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identity Eq. i|12|) implies that 
-^,(k)(^^ + ^) =A iifc (.,k). (17) 

There are some general consequences one can draw from 
this identity. First of all, in models of stationary forced 
turbulence Eq. J[J , the effective viscosity can not depend 
on time (nor will its FT depend on the frequency lo), and 
thus at most the renormalized viscosity could depend on 
wavenumber: v — v(k 2 ). The same holds for the "mass" 
term £ = £(fc), so that for stationary random forcing 
Eq. (|17|l immediately implies that Ajjfc(aj, k) = 0, and 

hence the vertex function in Eq. i|16|) is not al- 

tered by perturbative corrections to the equation of mo- 
tion. When Aijk(u), k) = 0, this vertex is that which 
corresponds to the convective term in the NSE Eq. Q). 
There is one important caveat here: this conclusion holds 
only for the case of zero-momentum transfer in the ver- 
tex (see Eq. dj). We are not claiming that the vertex 
cannot have nontrivial renormalization for finite momen- 
tum transfer. This result of vanishing vertex correction 
Ajjfc(w, k) is valid only for large-scale flows and is inde- 
pendent of the value of Ki. 

In conclusion, when we come to reconcile Galilean 
invariance with fluctuation phenomena, we know that 
the enforcement of this symmetry does impose nontriv- 
ial conditions on some of the physical parameters ap- 
pearing in the equation of motion. For example, re- 
cursive renormalization group (RG) procedures applied 
to the NSE explicitly generate higher-order nonlineari- 
ties in the renormalized fluid equation of motion. These 
higher-order nonlinearities, such as triple velo city prod- 
ucts, do however preserve Galilean invariance 

La- Cor- 
relation functions provide further information that can 
not be ascertained at the level of the evolution equation 
of motion alone. Here it is worthwhile to recall that the 
original operator-based perturbation theory developed by 
MSR requires three different vertices. They claim that 
the whole problem of strong turbulence is contained in 
a proper treatment of the vertex renormalization |l2j . 
The two additional vertices introduced by MSR have no 
counterpart at the level of the NSE. 

In contrast to the NSE, correlation functions can not 



be neatly separated into mean and fluctuating parts. 
This in general means constraints from Galilean invari- 
ance for the mean velocity component also can have im- 
plications for the fluctuating component. In particular, 
this leads to the WT identity Eq. 1|17[) . which provides 
a nontrivial relation between the fluid response function 
and the nonlinear vertex term (convective term). This 
dependence holds of course only for the large scale flow 
regime and says nothing about the small scale flow be- 
havior i.e., the inertial and dissipation ranges). This is 
so because (|17f) follows directly from l|12|l . which as can 
be checked, involves the vertex only for zero momentum 
transfer. 

In , the claim is made that Galilean invariance pro- 
vides no constraint whatsoever on the vertex renormal- 
ization, in sharp contrast to the results in These 
previous papers that examined the WT identities did so 
in the frame at zero mean velocity in which it would be 
impossible to test the claims of p}, which depend cru- 
cially on having Reynolds decomposed the fluid velocity. 
To test these claims requires working in a frame of gen- 
eral nonzero mean velocity. This paper has done this 
and demonstrated that although the claims in [7j have 
some validity, the fluctuating velocity component is still 
partly constrained, in particular in the limit of zero mo- 
mentum transfer. Thus our result provides a resolution 
to the seemingly differing conclusions drawn in Q and in 
0EJE3| In the process of addressing this issue, our paper 
has also provided the most general form of the WT iden- 
tities for the randomly forced NSE. These results can not 
be derived from the other treatments of the WT identities 
m 0- S @ by simply performing a Galilean transforma- 
tion of their results. Further understanding about the 
implications of Galilean invariance in the NSE can per- 
haps be gained by examining the analogy between this 
symmetry and global gauge invariance in QFT and ex- 
ploring the implications of gauge fixing procedures. 
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